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Abstract 

We derive the uncertainty relation for a quantum open system comprised of a Brow- 
nian particle interacting with a bath of quantum oscillators at finite temperature. We 
examine how the quantum and thermal fluctuations of the environment contribute to the 
uncertainty in the canonical variables of the system. We show that upon contact with the 
bath (assumed ohmic in this paper) the system evolves from a quantum-dominated state to 
a thermal-dominated state in a time which is the same as the decoherence time in similar 
models in the discussion of quantum to classical transition. This offers some insight into 
the physical mechanisms involved in the environment-induced decoherence process. We 
obtain closed analytic expressions for this generalized uncertainty relation under the con- 
ditions of high temperature and weak damping separately. We also consider under these 
conditions an arbitrarily-squeezed initial state and show how the squeeze parameter enters 
in the generalized uncertainty relation. Using these results we examine the transition of 
the system from a quantum pure state to a nonequilibrium quantum statistical state and 
to an equilibrium quantum statistical state. The three stages are marked by the decoher- 
ence time and the relaxation time respectively. With these observations we explicate the 
physical conditions when the two basic postulates of quantum statistical mechanics become 
valid. We also comment on the inappropriateness in the usage of the word classicality in 
many decoherence studies of quantum to classical transition. 
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Introduction 

This paper aims at expounding the relation of quantum and thermal fluctuations and 
their relative importance in aflFecting the behavior of a quantum open system [1]. The 
demarkation of quantum, classical and thermal regimes is not always clearly noted, their 
usual definitions or usage oftentimes are imprecise and misleading. For example, one hears 
the vague identification of a high temperature regime as the classical regime. One also hears 
the loose statement that classical is defined as the regime when the Heisenberg uncertainty 
principle ceases to hold. Many recent discussions of decoherence via interaction with an 
environment [2] view the disappearance of the off-diagonal components of a reduced density 
matrix in some special basis or the selection or emergence of a set of classical histories [3] 
as signalling a quantum to classical transition. What are the roles played by quantum and 
thermal fiuctuations in all these processes and how are these issues related to each other? 
[4] 

To seek an answer to these questions we focus on the derivation of a generalized 
uncertainty principle valid at finite temperature. The uncertainty in the canonical variables 
are caused by both the vacuum and thermal flucatuations of the system and the bath. We 
study this problem because while simple enough to yield analytic results it also embodies 
these issues fully. A summary of the main results of this investigation was recorded earlier 
in [5]. 

It is a well-known fact in quantum mechanics that a lower bound exists in the product 
of the variances of pairs of noncommutative observables. Taking the coordinate x and 
momentum p as examples, the Heisenberg uncertainty principle states that with (Aa;)^ = < 
> — < X the uncertainty function is 

Uq^ — (Ax)^(Ap)^ > — {T — 0, quantum mechanics) (1-1) 

The existence of quantum fiuctuations is a verified basic physical phenomenon. The origin 
of the uncertainty relation can be attributed as a mathematical property of Fourier analysis 
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[6] which describes quantum mechanics as a wave theory. Recent years have seen effort 
in estabhshing a stronger relation based on information-theoretical considerations [7]. In 
particular we refer to the papers of Anderson and Halliwell [8] and Halliwell [9] . 

In realistic conditions quantum systems are often prepared and studied at finite tem- 
peratures where thermal fluctuations permeate. At high temperatures the equipartition 
theorem of classical statistical mechanics imparts for each degree of freedom an uncer- 
tainty of kT/2. Thus the uncertainty function for a one-dimensional harmonic oscillator 
approaches the limit 

kT 

~ {high T, classical statistical mechanics) (1-2) 

where Q is its natural frequency Q. This result, obtained by assuming that the system 
obeys the Maxwell-Boltzmann distribution, is usually regarded as the classical limit. For 
a system of bosons in equilibrium at temperature T, the application of canonical ensemble 
gives the result in quantum statistical mechanics as 

h'^ hVt 

U^^ = — r[coth(— — )]^ (all T, quantum statistical mechanics) (1-3) 
4 2kT 

which interpolates between the two results (1.1) and (1.2) at T = and T » hVl/k. This 
result applies to a system already in equilibrium at temperature T. 

Under equilibrium and stationary conditions the finite-temperature uncertainty prin- 
ciple is just this simple. It has been studied before via thermofield dynamics [10]. Here 
we aim to tackle the fully nonequilibrium problem. Using the microdynamics of a quan- 
tum system as starting point we view thermal fluctuations as statistical variations of the 
coarse-grained environmental variables with which the quantum system interacts, the ex- 
act microdynamics of the system and the environment obeying only the laws of quantum 
mechanics. 

Our model is that of a collection of coupled harmonic oscillators where one is dis- 
tinguished as the system of interest and the rest as bath. This model has been studied 
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extensively before [2] [11-17]. We use the influence functional method [12] to incorporate 
the statistical effect of the bath on the system. At time to we put the system in contact 
with a heat bath at temperature T and follow its time evolution. We consider an initial 
Gaussian wave packet and calculate the spread < Ax < Ap >^ and the uncertainty 
function Ut due to both quantum and thermal fluctuations. The main result is given for- 
mally by Eq. (2.3.20), which, under the conditions of high temperature and weak damping, 
simplifles to Eq. (3.3.4) and Eq. (3.4.6) respectively. The results for an initial coherent 
or minimum-uncertainty state was flrst reported in [5]. Anderson and Halliwell [8] have 
recently shown that this result is equal to the lower bound of an information-theoretical 
definition of uncertainty with different initial coherent states. Here we derive the uncer- 
tainty relation for the general case of an arbitrarily-squeezed initial state where the squeeze 
parameter appears explicitly [17]. 

The paper is organized as follows: In Sec. 2, we begin with a short summary of the 
influence functional formalism applied to the Brownian motion model. (Details can be 
found in [12-15]. Readers familiar with this can skip to Sec. 2.2.) We then derive the 
reduced density matrix for a Gaussian wavepacket with nonzero initial position and mo- 
mentum and arrive at a finite temperature uncertainty function (2.3.20) in terms of the 
elementary functions ttj and coefficients introduced in [15]. In Sec. 3 we assume an 
ohmic bath and derive the forms of certain coefficients fn, Qn related to Ui, aij, expressing 
them in terms of the squeeze parameter S. We then work out the uncertainty function for 
a closed system (zero-coupling with the environment) of squeezed states (3.2.5). In Sec 
3.3 and 3.4 we obtain the uncertainty function for high temperature and weak damping 
cases, given respectively by (3.3.4) and (3.4.6). For weak couplings, if the initial state is a 
minimum-uncertainty state, we recover the result (3.4.7) obtained earlier by us [5] and by 
Anderson and Halliwell [8] . In Sec. 4 we use our results to discuss the relative importance 
of quantum and thermal fiuctuations and show how this bears on two related issues: One 
on the relation of uncertainty to decoherence studies [2] , and the other on the demarkation 
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of quantum, classical and thermal regimes. In our results one can identify a time where 
the thermal fluctuation overtakes the quantum fluctuation. This time is identical to the 
decoherence time obtained in earlier studies of quantum to classical transition [2]. The 
other time scale of significance is that of relaxation time. We caution that after the de- 
coherence time, although the system is describable in terms of probabilities, it cannot yet 
be regarded as classical because of the spin-statistics effect, and it remains in a state of 
nonequilibrium. After the relaxation time the two basic postulates of statistical mechanics 
[18] become valid and one can depict the system under the premises of equilibrium (but 
still quantum) statistical mechanics. The system has to be in a sufficiently high temper- 
ature when the spin-statistics can be adequately represented by the Maxwell- Boltzmann 
distribution (particle becomes distinguishable) for it to be aptly called classical. 
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2. Quantum and Thermal Fluctuations in an Open System 



2.1 Influence Functional for a Brownian Particle in a Harmonic Oscillator Bath 

Consider a Brownian particle with mass M and natural frequency Qq. Its environment 
is modeled by a set of harmonic oscillators with mass m„ and natural frequency a;„. The 
Brownian particle is coupled linearly to the nth oscillator with strength Cn- The action of 
the combined system and environment is 



S[x, q] = S[x] + Sb[q] + Sint[x, q] 
t 

= I ds 



Q L n n _ 

(2.1.1) 

where x and qn are the coordinates of the Brownian particle and the n-th bath oscillators 
respectively. We are interested in how the environment affects the system in some aver- 
aged way. The quantity containing this information is the reduced density matrix of the 
system Pr{x, x' , t) obtained from the full density operator of the system and environment 
p{x, q; x', q'; t) by tracing out the environmental degrees of freedom (g„, g^) 

+00 +00 

Pr{x,x',t) J dqn J dq'^ p{x,q;x',q',t)S{qn-q'J (2.1.2) 

" -oo -oo 

The reduced density matrix evolves under the action of the evolutionary operator Jr in 
the following way: 

+00 +00 

Pr{x,x',t)— J dxi J dx'^ Jr{x,x' ,t \ Xi,x'j^,0) Pr{Xi,x'^,0 ) (2.1.3) 
— oo — oo 

In general, this is a very complicated expression since Jr depends on the initial state. If 
we assume that at a given time t = the system and the environment are uncorrelated, 
i.e. that 



(2.1.4) 



then the evolution operator for the reduced density matrix can be written as 



Jr{xf,x'f,t I Xi,Xi,0) ^ J Dx J Dx' exp^|5[a;] F[x,x'] (2.1.5) 

Xi x'. 

where F[x,x'] is the Feynman- Vernon influence functional [12]. If the environment is 
initially in thermal equilibrium at temperature T, for the problem described by (2.1.1), 
the influence functional can be computed exactly. The result is well known [12-15]: 



t Sl 



F[a;, a;'] = exp< ~ ^ / '^"^i / ^"^2 x{si) - x'{si)^r]{si - S2) x{s2) + x'{s2)^ 





t Sl 



(2.1.6) 



n 



dsi / ds2 



x{si)-x'{si) iy{si-S2) x{s2)-x'{s2) 





The non-local kernels r]{s) and ^{s) are deflned as 



and 



where 



+00 

, , [ du ^, . ^ Tiuj 
^\^) — I — /(a;) coth — — cosa;s 
J 71 2kT 



r^is) = ^ 7(^) 
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+ CXD 



duj I{co) 



cos LOS 



(2.1.7a) 



(2.1.76) 



(2.1.7c) 



Here /(a;) is the spectral density function of the environment, 



/(a;) = 5^5(u;-a;„); ""^^ 



2mna;r 



(2.1.8) 
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The real and imaginary parts of the exponent in the influence functional F[x, x'] are 
understood as responsible for dissipation and noise respectively, thus the names dissipation 
and noise kernels are given to 'r]{s) and The most general environment would thus 

engender nonlocal dissipation and colored noise. (We refer the reader to Ref [1, 15] for 
a discussion of the fluctuation-dissipation relation and the time scales of the relevant 
processes.) 

An example of the spectral density is given by [13] 

I{u;) = M^oouC^Ye-^ (2.1.9) 

where is a frequency scale usually taken to be the cut-off frequency A. The environment 
with this kind of spectral density is classified as 



ohmic if s = 1 
supra — ohmic if s > 1 
suhohmic if s <1 



(2.1.10a) 
(2.1.106) 
(2.1.10c) 



The propagator (2.1.5) has been calculated before [12-15]. (We use the notation of 
Ref. 15) 



Jr{xf,x'f,t\ Xi,xl,0) = Zo{t)exp-^ ui{0)Xi + U2{0)Xf 



Here we have introduced 



+ i 



ui{t)Xi + U2{t)Xf Yf + iaii{t)Yi 



ai2{t) + a2i{t) YiYf + ia22{t)Yf 



Xi = ^{xi + x'i) 



(2.1.11) 



(2.1.12a) 



and 



2^"-' ' (2.1.126) 
Yf = x'f- Xf 

The elementary functions Ui{s),i = 1, 2 in (2.1.11) satisify the following integro-differential 



equation 

with boundary conditions 



s 

+ 2 J ds'r]{{s - s')ui{s') + nlui{s) = (2.1.13) 





„.(0)=i U.(0)=o 

«i(0) = I «2(0) =1 ^ ^ 



The coefficient functions Oy(t), i,j — 1,2 are defined by 



t Si 



a,At) = jds,jds,v,M.(s,-s,)vAs.) (2.1.15) 





where 



vi{s) = U2{t - s) (2.1.16a) 

and 

V2{s) = ux{t - s) (2.1.166) 
It is clear that the functions Ui{t) and aij{t) only depend on the kernels 'q{s) and ^{s). 
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2.2 Reduced Density Matrix For A Gaussian Wavepacket 

We now consider a Brownian harmonic oscillator with an initial wave function (Calde- 
ria and Leggett 1985, Unruh and Zurek 1987 in [2], Hu, Paz and Zhang 1992 in [15], Paz, 
Habib and Zurek 1992 in [16]) 

0) = ^/N~o exp{-^^^^ + ^pox} (2.2.1) 

where a is the initial spread, and xq and po are the averaged initial position and momentum 
of the Gaussian wave packet. The initial reduced density matrix is given by 



Pr{xi, a^j, 0) = No exp< -poa^i + ^Po^i > 



(2.2.2) 



where 



No = Noexp[-^] (2.2.3) 
Substituting this and the propogator (2.1.11) into (2.1.3), one would get 



+ 00 +00 

Pr{xf,x'f,t) ^ ZoNoexp^^^^-U2{t)XfYf + ia22{t)Y}^ X dX dY 



X exp< 



— oo —oo 

l(X,Y (H^^{t) H^2{t)\(xA , (h^{t)Y (X, 



2\YJ \H2i{t) H22{t) J \YJ ' \b2{t) J \Yi 

(2.2.4) 



where 



h2{t) = ^U2{Q)Xf - ^ [ai2(t) + a2i{t)]Yf + ^po 



(2.2.5a) 
(2.2.56) 
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and the matrix elements of H(t) are given by 





1 




(2.2.6)a 


Hi2{t) = 


H2l{t) 




(2.2.66) 


H22{t) = 


1 

4^2 ^ 


j^a22{t) 


(2.2.6c) 



Performing the Gaussian integrals over X and Y in (2.2.4), one gets 



= \ {-A itHi ^ Hi)- ^"'<"^^ 



where 



^'^^^ = ndetH(t) {~'^^^^^^ + ^''^^^ t""''^*^ + I''" 
" a2ndetH(t) {"^^^^^^^^^^ + '^''^^^ t""''^^^ + 



The prefactor 



2\Yf J \Q2i{t) Q22{t) J \Yf J ' \B2{t) J \Yf 

(2.2.7) 



= -;tf^t\^o - ^#^.0, (2.2.8a) 
^ n(72detH(t) " n^detH(t)^ ' ^ ' 



(2.2.86) 



(2.2.9) 

depends only on the parameters xq, Po and time. The matrix elements of Q(t) are given 

by 

gn(0 = (2.2.10a) 
nV2detH(t) ^ ^ 

12 



h 



r detH(t) [a 



ai2(t) + 021 (t) 



- ^wi(0)tti(t) } (2.2.106) 



Q22(t) = T«22(t) + 



h 

2 



det H(t) 



1 r i2 
ai2(t) + a2i(t)|wi(0)?ii(t) - ai2{t) + a^iit) 



(2.2.10c) 



2.3 Finite Temperature Quantum Uncertainty Function 

To calculate the averages of observables, it is convenient to use the Wigner function 
associated with the reduced density matrix pr defined as 



+ 00 



w. 



{x,p,t)= j dY e*^'^p,(x-^y,x+^y,t), 



with an inverse relation given by 



p^{X--Y,X+-Y,t) 



1 f dp 



(2.3.1a) 



+ 00 

J ^^'^""^ ^r{X,p,t) (2.3.16) 



Applying the above Wigner transform to the reduced density matrix (2.2.7), we get 



+00 

{X,p,t)= J dY e*P^ Zo{t) 



^"^P^ 2[Yf) [Q,,{t) Q22{t) ) {Yf ) + [B,{t) ) [Yf 

(2.3.2) 

where the superscript T denotes the transpose of the matrix. The quantum statistical 
averages of an observable of the system, e.g., x", or with respect to the reduced density 
matrix Pr{t) are given by 
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+ 00 +00 



<X^>T^ I dx I ^X^Wr{X,p,t) 



27rh 



(2.3.3) 



-oo — oo 



and 



+00 +00 



<p">T= / dx I ^p-Wr{X,p,t) 



(2.3.4) 



—00 —00 

where the subscript T indicates that the environment is at temperature T. Obviously the 
averages have both quantum and thermal contributions. It is easy to show that 



and 



thus 



Similarly, 



and 



thus 



< X >T-- 



Bijt) 
Quit) 



(2.3.5) 



<X>T 



(2.3.6) 



iAxr = <x^>r-<x>^,= ^ 



(2.3.7) 



detQ(t)] 5 

Qii{t) 



B2{t)-Bi{t) 



Ql2{t) 

Quit) 



(2.3.8) 



2 ^ 2 det Q 



< 



P >t] 



(2.3.9) 



(Ap)| 



< p^ >T - < P >T 



detQ 
Qii(t) 



(2.3.10) 
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The finite temperature quantum uncertainty function defined as 



(2.3.11) 



follows: 



= h 



2 dctQ(t) 

2 1 Q22{t) 



Ql2(t) 



Using (2.2.10a-c), one can rewrite 



Ql2(t) 

Quit) 



= ifo{t)+ifi{t) 



and 



(2.3.12) 



(2.3.13) 



Q22{t) 



9o{t) + gi{t) + g2{t) 



(2.3.14) 



Quit) 

where the functions fnidn denote terms of order n in temperature T. Explicitly they are 
given by: 



Wl(0)w2(t) -Wl(t)w2(0) 



(2.3.15) 



and 



hit) 



^ ''2^^an(t) + ai2(t) + a2i(t) 



«2(0)i tt2(0) 



(2.3.16) 



9oit) = 



uiit) 

«2(0) 



(2.3.17) 



2^T^ I 

= 77—7772 1 Vi(t) +wi(0)wi(t) 



X 



n[^2(0)]' 

ai2it) + a2i(t) 



+ 



+[«i(0)]' a22(t) , 



(2.3.18) 
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92{t) = 



[MO)]- 



4:aii{t)a22{t) - [ai2{t) + a2i{t)\ 



(2.3.19) 



The finite temperature uncertainty function is then formally given by 



UT{t)= go{t) + gi{t) + g2{t) + Mt) + h{t) 



(2.3.20) 
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3. Uncertainty Principle and Squeezed States 



3.1 Finite Temperature Ohmic Environment 

We now proceed to analyze the finite temperature uncertainty function for an ohmic 
environment. In the cases of high temperature and weak couplings, one can obtain rela- 
tively simple analytic expressions for the uncertainty relation for a general squeezed states. 
Results for the (unsqueezed) coherent state or the minimal uncertainty state were obtained 
by Hu and Zhang [5] earlier. For simplicity, we set the mass of the Brownian harmonic 
oscillator M = 1. The spectral density for ohmic dissipation is 



2 



/(a;) = 7oa; e 



(3.1.1) 



With this, the dissipation kernel becomes 




(3.1.2) 



If we assume that the cutoff frequency A is very large, then 




7(s) ~ 7o(^(s) 



(3.1.3) 



so the damping becomes local. However, the noise kernel 




(3.1.4) 







is generally nonlocal except at very high temperatures. 




(3.1.5) 
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whence 



v{s) ~ -j^io5{s) (3.1.6) 



corresponding to white noise. 

For local dissipation, equation (2.1.13) becomes 



d^Ui{s) , dui{s) , 
+ 7o — J 1- 



nl - 27o5(0) Ui{s) = -2S{s)ui{0) (3.1.7) 



ds^ ds 

The boundary term on the right side of equation (3.1.7) can be neglected since it gives no 
contribution to the solutions with the boundary conditions (2.1.14). We also can redefine 
the natural frequency by 

{j]2_27o5(0)}^J]2^ (3.1.8) 

then equation (3.1.7) becomes the familiar equation of motion for a damped harmonic 
oscillator. Let us define a damping parameter 



and the effective frequency [15] 



= (^^0 - (3.1.10) 
Then in the case of underdamping, a < 1 the elementary functions simplify to 



= sinn(t-^ 

^ ^ sinllt ^ ^ 

and 

u^^s) = e^To(^-^) (3.1.116) 
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^ COS Qt 1 

—si 7o 

sin lit 2 



sin Qt 



hot 



The derivatives of these elementary functions are 

«i(0) 

and 

^2(0) ■ 

smOt 

• / \ cosOt 1 
U2{t) = ii ■ ^ , - -7o 
smiZot 2 

By using Eq. (3.1.12a-d), it can be shown that 



(3.1.12a) 
(3.1.126) 

(3.1.12c) 
(3.1. 12d) 



^ — ^ sin 2nt-a ^-i^ sin^ nt \ 6"^°* 
4:5 25 J 



1 + 52 



and 



(3.1.13) 



9oit) = ^e-^^°* 
where we have introduced a squeeze parameter 



(3.1.14) 



which measures the spread in the initial Gaussian wavepacket. {5 
coherent state or minimal-uncertainty state.) 
It can be shown that (see Appendix) 



(3.1.15) 

1 corresponds to a 



(3.1.16) 



9i{t) 



+ 



+ 



.Oo r 

^2^[ 
5 



2Ql 



4nl 



li_7oO 
^5^1 Ql 

1-1-2^" 

52 2np 



sin 2Qt — 



1 - ^ - 



52 2Ql 



cos 2Qt > [ss 



sin 20,t — cos 20,t )■ \sc] 



1 + 



1 

51 



+ ^sin2nt + 



1 - 1 - ^ 

52 21^2 J 



COS 2Q,t > \cc 



(3.1.17) 
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and 



92{t) 




(3.1.18) 



where the functions [ss], [sc] and [cc] are defined in (A.l) 
3.2 Zero Coupling Limit 

Let us first examine the simplest case of zero-couphng. It corresponds to an isolated 
harmonic oscillator taken as a closed quantum system. We expect to recover the familiar 
results in quantum mechanics. 

Assuming 70 = 0, then 7(5) = and i'{s) = 0. Therefore, 



f,{s)=g,{t)=g2{t) = 



(3.2.1) 



and 



n = a 







(3.2.2) 



From Eq. (3.1.14) and (3.1.15) one gets 



fo{t) 



4S 



sin 2^0^ 



(3.2.3) 



and 



go{t) 



1 



(3.2.4) 
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We find the quantum uncertainty function for an initial squeezed states to be 




(3.2.5) 
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The time-dependent term is the result of quantum dispersion, 
herent state, 

S=l 



For the (unsqueezed) co- 
(3.2.6) 



we recover the Heisenberg uncertainty relation 

Urit) = ^ (3.2.7) 
With this we can also understand why 5 = 1 is called a minimun-uncertainty state. 

3.3 High Temperature Limit 

As one can see from (3.1.3) and (3.1.6), at high temperatures both the noise kernel 
and the dissipation kernel for ohmic dissipation become local. In this limit the functions 
/ and g simplify to (See Appendix) 

fi(t) = aro sin^ nte-^°* (3.3.1) 



5'l(^) = To' 





1 - e-^°^ 


\ 45 





+ ^^o? sin^ 9.t —OL sin 2^t l>e"^°* (3.3.2a) 



25 



45 



and 



92{t)^Tl{\ 



\^o? sin^ 



(3.3.26) 



where we have introduced the dimensionless parameter 

TO = = - 3.3.3 

We will also use e{uS) and e to denote the same quantities as eo, but with a; and O replacing 
r^o in (3.3.3) respectively. The uncertainty function for ohmic dissipation at the high 
temperature limit is then given by 
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n 



1 - e-To* 



+ T0 



l-e 



-lot 



-lot 



- ToQ;< 


f 1-^2 


sin 2f]t + Q! 


To - 


1 + 52- 


sin^ l]t j 




25 . 






sin 2Vlt + a 


To - 


1 + 5^1 


sin^ l^tj 


45 


25 \ 



(3.3.4) 



This is the first main result of this paper. 

Alternatively, since the noise kernel is also a Dirac delta function, from (2.1.15), we 

have 



kT 

aij{t) ~ —70 



t 

J ds Vi{s)vj 



is) 



(3.3.5) 



They are easily computed to be 



aii(t) ~ 



kT 



3^"*|l- ^ -a^cos2nt + asin2nty-'^°*^ (3.3.6a) 



cos — q; sin Vtt 



+ 



cos + a sin Vtt 



-lot 



(3.3.66) 



022 (^) 



kT 



q2 

- CK^ cos 2nt - a sin 2^^ 



-lot 



(3.3.6c) 



Substituting these and (3.1.12a-d) into (2.3.16) and (2.3.18-20), we can also arrive at 
Eq. (3.3.2-4). 

We see that there are two factors at play here: time and temperature. Time is 
measured in units of the relaxation time proportional to trei = lo^i and temperature is 
measured with reference to the ground state energy %Q,q/2 of the system. There are also 
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two parameters involved, a and 8. Let us now take a closer look at this high temperature 
uncertainty function under different conditions: 

a) At t = 0, we find 



t/T(0) = ^ (3.3.7) 



At t » 7o \ we find 



UHt) ^ (-) (3.3.8) 

These comply with the expected results from quantum mechanics and classical statistical 
mechanics as stated in the Introduction. 

b) At short times t « ^q^, we can expand the uncertainty function (3.3.4) as 

C/y(t) ~ _[l + 2(ro5-l)7ot + C'(t')] (3.3.9) 

Prom this we can identify the time when thermal fluctuations overtake quantum fluctua- 
tions, i.e., 

ti ~ — — — r (high temperature) (3.3.10) 
4:kT^oo 

which is seen to depend on S, the squeeze parameter of the initial state. This time is 
identical to the decoherence time derived in [2] [15, 16]. That the transition time from 
quantum to thermal dominance in the uncertainty relation is related to the decoherence 
time scale was first noted by the present authors in [5] . 

c) For a minimum uncertainty initial state (S = 1), 
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6--^°* + To 



1 - e-^°* 



Toa^(To - l)sin2 



d) For weak damping, one can neglect all 70/^^0 terms and find 



-lot 



+ ro 
+ 



45 
16^2 



-Jot 



--lot 



sin^ 9.te-^^°* 



e) For weak damping, and for an initial minimum-uncertainty state. 



■Jot 



1 - e-^°* 



(3.3.11) 



(3.3.12) 



(3.3.13) 



It is seen that the first term is of purely quantum nature whereas the second term is of 
thermal nature, their contributions to the uncertainty of the system arise from quantum 
and thermal fluctuations respectively. 



3.4 Weak Damping Limit 

At arbitrary but finite temperatures, the noise kernel (3.1.4) is not a Dirac delta 
function, and in general there is no closed form for the uncertainty function. However, if 
the damping is weak, we can find an approximate expression for the uncertainty function. 

Assuming that 

a^^«l (3.4.1) 

then 



1 „ 
2 



(i7o)2 + (f)-a;)2 



(3.4.2) 
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which suggests that the major contributions to the integrals in (A.4a-c) come from a small 
region near u = il. Using Eqs. (A. 14) for (3.1.16-18), we obtain 



flit) 



70 

2no 



coth er- sin^ Qt e"'^°* 



9i{t) -cothe^ 



1 + 
26 



1-e 



-lot 



25 



a sin 2nt \ e"^°* 



and 



g2{t)^ [^cothe, (l-e-2To*)l 



(3.4.3) 
(3.4.4) 

(3.4.5) 



The uncertainty function (2.3.20) is then given by 



1 1 

■^UT(t)^ - { e-^o* + cothe. 



-lot 



+ coth e 



(1-5) 
A5 



2 r 



1 



+ <! — -— sm 2Qt + a 



-lot 



coth er 



46 
25 



a sm.2Q,t je 
sin^ Qt\- e 



-lot 



-270* 



(3.4.6) 



This is the second main result of this paper. One can deduce various limits from this 
expression. For example, it is obvious that again at t = 0, when the initial uncorrelated 
conditions is assumed valid, Ut{0) = which is the Heisenberg relation (1.1). At very 

long time (t » ^q^), Uxit) is insensitive to 5 and approaches U^^ as in (1.3) at finite 
temperature (as VL ~ ^q). That means the Brownian particle approaches an equilibrium 
quantum statistical system. (For supraohmic bath this may not always be true). We can 
also see that for a T = bath (cothe = 1), Urit) has a leading term given by /A 
(the Heisenberg relation) followed, for squeezed states 5 ^ 1, by terms of order and 
a depicting both decay and oscillatory behavior. This is, of course, due to the action of 
quantum fluctuations alone. 

For a minimum- uncertainty initial state {5 — 1), we get for all finite temperatures 

-2 r _ _ ^ 2 



e"^°* + cothe 



1 - e-^°* 



+0{a^) {weak coupling) 



(3.4.7) 
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This result was obtained in [5] earlier. Obviously, taking the high temperature limit of 
(3.4.7) reduces to (3.3.13) above. 

Notice that in (3.4.7) there is no linear order damping term. In this calculation we 
introduced a cutoff frequency A in the spectral density /(a;) ~ a; exp(— a;^/A^) which 
leads to a divergent term dependent on A. It is removed by a frequency renormalization 
procedure standard in field theory. (See Eq. (A. 12c) and discussions in the 1992 paper of 
[15] and Appendix B of the 1983 paper of [13]). 

At short times (t « ^q^), 



2 

Urit) ^ ^ [l + 2(5 coth ^ - l)7ot + 0{t')] (3.4.8) 

This simple expression is revealing in several aspects: The first term is the ubiquitous 
quantum fluctuation, the second term is the thermal contribution, which depends on the 
initial spread and increases with increasing dissipation and temperature. The time when 
thermal fluctuations overtake quantum fluctuations is (assuming the temperature is higher 
than the ground state energy): 

= 3—- (3.4.9) 

27o(5coth^-l) ^ ^ 

Even though this result has not appeared in decoherence studies of quantum to classical 
transitions [15], we expect this to be equal to the decoherence time scale tdec calculated 
for weak coupling for all temperatures. 
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4. Discussion 

Combining the generalized uncertainty relations obtained here with the recent findings 
of the environment-induced decoherence studies based on the same open-system frame- 
work, where two characteristic times — the decoherence time tdec and the relaxation time 
trei — are defined, one can address two basic issues in quantum mechanics and in (nonequi- 
librium) quantum statistical mechanics i.e., 1) expound the relation between quantum, 
thermal and classical: more specifically, depict the role of quantum and thermal fluctua- 
tions in the quantum to classical transition; and 2) explicate the physical conditions for 
the realization of the two basic tenets of (equilibrium) quantum statistical mechanics from 
quantum dynamics. 

Quantum statistical mechanics of a macroscopic system is derived from the quantum 
dynamics of its microscopic constituents under two basic postulates [18]: i) random phase, 
and ii) equal a priori probability. The first condition enables one to assign probability 
distributions to a system occupying certain quantum states. It requires the suppression of 
interference terms in the wave function or that the reduced density matrix of the system 
be approximately diagonal. The second condition is the basis for the microcanonical 
ensemble. When applied to the stipulations of a canonical ensemble, it ensures that under 
general conditions the system when put in contact with a large bath would eventually 
equilibrate with it. We want to examine the processes by which these two conditions are 
attained from a more basic level, starting with the microdynamics of a system of quantum 
particles. Specifically, we want to see if there is a characteristic time when the phase 
information is lost (Postulate i) and another time when the system attains equilibrium 
with its surrounding so that all accesible states (in the combined system and environment) 
are equally probable (Postulate ii). 

Now for the first issue which embodies the question we raised in the beginning and 
which underlies the premises of the second issue: On the one hand, one often regards the 
regime when thermal fluctuations begin to surpass quantum fluctuations as the transition 
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point from quantum to classical. On the other hand, from the wave picture of quantum 
mechanics we know that a necessary condition for a system to behave classically is that 
the interference terms in its wave function have to diminish below a certain level, so that 
probability can be assigned to classical events [2] or that classical decoherent histories can 
be well-defined [3]. This is known as the decoherence process. One should ask if there 
exists any relation between these two criteria of classicality - one based on the uncertainty 
relation and the other on decoherence. In particular we want to see how the uncertainty 
function Urit) changes from the initial quantum fluctuation-dominated condition to a 
later thermal fluctuation-dominated condition. In Eq. (3.3.13) we show that under the 
stipulated conditions they are indeed equivalent: The time the quantum system 'decoheres' 
is also the time when thermal fluctuation overtakes quantum fluctuations. 

However, we issue a warning here that the regime after thermal fluctuations dominate 
should not be called classical, as is customary in many quantum to classical transition 
studies. In fact, after the decoherence time only the first postulate of quantum statisti- 
cal mechanics (QSM) is satisfied, the system can be described by non- equilibrium QSM. 
Only after the relaxation time, when the second postulate is satisfied, can the system 
be correctly described by equilibrium QSM. The classical regime starts at a much later 
time and it would have to be at a sufficiently high temperature. It is well-known that 
quantum statistical eff'ects can be important at very high temperatures (e.g., Fermi tem- 
perature for metals) . This is due to exchange interactions of identical particles, a distinctly 
quantum effect. Only when the statistical properties of fermions and bosons can be approx- 
imated by distinguishable particles, usually at high temperatures when the Fermi-Dirac or 
Bose-Einstein statistics approaches the Maxwell-Boltzmann statistics, can the system be 
rightfully called classical. In this regard quantum carries two meanings, one refers to the 
interference effect and the other refers to spin-statistics effect. 

The analysis of nonohmic environments will be reported later [19]. There one expects 
to see some qualitatively new behavior in the case of supraohmic spectral density function 
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and at low temperatures, when quantum fluctuations can play a more determinant role 
which can make the decoherence time longer and the relaxation of the system to equilibrium 
more difiicult to achieve. 
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Appendix 



Let us define 



then 



t t 



[ss] = ^ /^^i/ smQsi smQs2 e-^^°^''+'^^ 





+ CXD 



X 7o / — a;cothe(a;) cosa;(si — S2) 

J TT 



t t 



[sc] = ^jdsij ds2 sinf]si cosl]s2 e-^^°^''+'^^ 



+00 

X 7o / — a; cothe(a;) cosa;(si — S2) 

J TT 




t t 

[cc] = ^JdsiJ ds2 cosflsi cosfisa e"5To(*i+*2) 



+00 

X 7o / — a;cothe(a;) cosa;(si — S2) 




(Ala) 



(A.lb) 



{A.lc) 



aii(t) = 



ai2(^) + a2i{t) = 2 — J — e^"'"^\smQt [sc] - cosQt [ss]\ 
sin Qt J 

0^22 (i) — — 5 — |cos^ Qt [ss] — sm2Qt [sc] + sin^ Qt [cc]\ 
sin Qt ^ J 



(yl.2a) 
(yl.2c) 



By switching the orders of the integrations in (A. la), (A.lb) and (A.lc), and by first 
performing the integrations of si and S2, we obtain 
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+ 00 

iooc[t) ^ - / — a;cothe(a;) 

Z I TT 



and 



ho 



U)± 



(a) In the high temperature limit 



+00 



f duo 




1 TT 


7- +7+ 



cos ijOt 





+ CXD 





+00 



kT 



^ , , kT f duj[ -1 . kT 



^ , , kT f duj[ -1 /cT 



In particular, 



/^e(0) ^ 



/c^c(O) ~ 



then 



\ss\ ~ 



1 - 



Q2 

- cos 2nt + a sin 2nt 
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r 1 kT 

[sc\ ~ — xT^as 1 - 



[cc] ~ 



kT 02 
X202 



cos 2Qt + a sin 20t 



-70 1 



1 + 



202 



Q2 , ^ 

^ + cos 20t - a sin 20tJ 6-"^°* \ 



{A.8b) 

(A.Sc) 



Substituting Eq. (A.8a-c) into Eq. (3.1.17), (3.1.18) and (3.3.19), we obtain Eq. (3.3.1) 
(3.3.2) and (3.3.3). 

(b) For weak damping, 

If the damping constant is very small, i.e.. 



a «1 



{A.9) 



then 



ho 



(AAO) 



which suggests that the major contributions to integrals (A. 4a), (A. 4b) and (A. 4c) come 
from a small region near u = Q. Thus we can make the following expansion: 



tanh tanh 



2kBT 



1 - 



1 - 



sinh er J O 



Q — OJ 



to get 



(All) 



Ijc{t) 



2 e"^^°* cosQt^ 1-a 



tanh e 



1 - 



sinh er 



sin Qt 
cos 



(^.12a) 



tanh e 



q; 



sinh er 



cos Ot 
sin fit 



tanh 



1 - 



sinh Cr 



The second term drops out when a high frequency cutoff A is assumed. 



lujcit) e-^To* sinOt<| 1 - a 

tanh er 



1 - 



sin Qt 



sinh er J cos Clt 



(yl.l26) 



(A.12c) 
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and 



/7c(0) 



1^ 



tanh er 



{A.lSa) 



_¥l_ 

tanh er 



a 



1 - 



sinh er - 



{A.13b) 



Substituting (A.12a-c) and (A.13a-b) into (A.3a-c), we obtain 



r 1 ^ n 1 

[ss\ ~ — cotner'^ 1 



1 + a sin 2nt^ 6"^°* 



(A.Ua) 



[sc] ~ — cotherCt-j 1 — 



cos 2flt + a sin 2Qt 



[cc] ~ — coth e^ ^ 1 — 



1 — q: sin 2Qt 



-Jot 



(A.14c) 
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